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Abstract. A model of a ladder of Josephson junctions in a magnetic field is 
considered. The topological features of the zero-temperature phase diagram, 
as a function of charging energy and small deviations from commensurabil- 
ity of the vortex lattice, are strongly dependent on the field. In addition to 
superconductor-insulator transitions, commensurate-incommensurate tran- 
sitions and depinning by quantum fluctuations of the vortex lattice are also 
possible. The critical behavior of the superconductor-insulator transition at 
/ = 1/2 is found to depend on the ratio between interchain and intrachain 
Josephson couplings and is in the universality class of the classical XY-Ising 
model. 



1. Introduction 

Arrays of Josephson junctions can currently be fabricated in any desired 
geometry and with well-controlled parameters. There has been consider- 
able interest in this system, both experimentally and theoretically, and 
as a tool to investigate several interesting effects like disorder, dissipation 
and quantum fluctuations [1]. At low temperatures where capacitive effects 
dominate, the array undergoes a superconductor to insulator transition as 
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a function of charging energy [2, 3, 4, 5, 6]. These charging effects arise 
from the small capacitance of the grains making up the array and leads 
to strong quantum fluctuations of the phase of the superconducting order 
parameter. The universality class of this T = superconductor-insulator 
transition is currently a problem of great interest specially in relation to 
experiments[2, 3] and theoretical predictions of universal properties[4, 7]. 
In fact, some of these properties, like the universal resistance, have been 
corroborated by measurents[3] and numerical calculations [8]. A Josephson- 
junction ladder provides the simplest one- dimensional version of an array 
in a magnetic and is an excellent system to study several interesting zero- 
temperature transitions due to commensurability effects of the flux lattice 
introduced by the field [14, 16]. 

A Josephson-junction ladder can be constructed by placing supercon- 
ducting islands at the sites of a one- dimensional ladder as indicated in Fig. 

1. At each site r of the ladder we associate a phase 9 r and charge 2en r , 
representing the superconducting element which is coupled to its neighbors 
by Josephson couplings E x or E y . The variables n r and 9 r are conjugate to 
each other. If only capacitance to the ground is considered the interaction 
Hamiltonian is given by the self-charging model 

r r <rr'> 

The first term in Eq. (1) describes quantum fluctuations induced by the 
charging energy E c = 4e 2 /C of a non-neutral superconducting grain, where 
e is the electronic charge and C the effective capacitance. The second term 
is the usual Josephson-junction coupling between nearest-neighbor grains. 
A rr t = (27r/<I> ) j r A - dr and A is the vector potential due to the external 
magnetic field and the gauge-invariant sum around a plaquette ^ A rr t = 
2irf with / = &/& . Here / is a dimensionless frustration parameter and 
<i> the flux quantum. Interesting quantum critical behavior can take place 
for rational values, / = p/q, of the frustration. Since the Hamiltonian is 
periodic under / — ► / + n (n integer), / = is equivalent to / = p/q with 
q=l. 

2. Commensurability effects 

As a result of the competition between the periodicity of the vortex lattice 
and the underlying pinning potential provided by the ladder, different phase 
transitions are possible as a function of the charging energy and the field 
[14]. As the magnetic field is increased from zero, a transition into a vortex 
state can occur where the magnetic flux first penetrates the ladder. This 
transition can be viewed as a commensurate-incommensurate transition 




Figure 1. Periodic Josephson-junction ladder. Superconducting islands (squares) are 
coupled to their nearest neighbors with inter-(_E x ) and interchain (E y ) coupling constants. 
$ is the magnetic flux through an elementary cell. 



described by the sine-Gordon Hamiltonian [16, 17] and is the ladder analog 
ol flux penetration in type-II superconductors. For small fields, the phases 
in different branches of the ladder are locked to each other while in the 
vortex state exponentially interacting kinks (vortices) appear that unlock 
the phases. Inclusion of charging effects leads to an additional insulating 
phase. In the vortex state, however, the vortex lattice can also become 
commensurate with the ladder at rational values of the flux quanta per 
plaquette f = p/q. The behavior of these commensurate vortex phases 
as a function of small deviations from commensur ability Sf = f — f and 
charging energy is an interesting problem to be addressed. One can expect 
that for / = p/q (p and q relative primes) the resulting phase diagram will 
have a strongly dependence on q. In fact, many different properties of an 
array in a field [4, 9] depend only on q. 

To study the global features of the phase diagram we use an effective 
action describing fluctuations from the commensurate state [14]. Starting 
from the path-integral representation of Eq. (1) and introducing auxiliary 
fields we obtain an effective action F of the form 

F = fdxf dT{±K x [(£<h + TrSf f + ( j-J 2 - xSf f] 

+ ±K T [(£<h) 2 + (£fo) 2 ] -wcosqfa - h)} (2) 

where cj)i(x), ^(x) are phases measuring deviations of the order param- 
eter in each branch of the ladder from the commensurate state, and 
K = K X K T ~ E y l E c and w ~ E x are effective couplings. 



4 



When 8f = 0, Eq. (2) is in the form of a Gaussian approximation of 
two coupled classical XY models which has been studied previously with E c 
playing the role of an effective temperature [fO, 11]. When vortices in (f>i, 4>2 
are included different critical behaviors result as a function of q. If q = 1, a 
single Kosterlitz-Thouless (KT) transition occurs as E c increases separating 
a commensurate (superconducting ) phase with long-range order in the 
phase difference <j)\ — 4>2 and algebraic order in (f>i,(f>2, from a disordered 
(insulating) phase where correlations decay exponentially. For q = 2, recent 
studies indicate that it is either nonuniversal or first order if it is a single 
transition and it is well described by the XY-Ising model [11, 18]. For q > 
an intermediate incommensurate (superconducting) phase is possible 
with algebraic order in all correlations. 

The effect of small deviations from commensurability, 8f > 0, can be 
studied by using a change of variables ip\ = 4>\ + 4>2, "02 = <f>i — 4>2 which 
leads to 

F = j dx j dri-lU^f + -iU^f 
J J A ox 4 or 

1 6 1 6 

+ -K x (—4> 2 - 2irSf) 2 + -K T {—i} 2 f - w cos # 2 } (3) 

and noting that this Hamiltonian has decoupled into a Gaussian in ipi and 
a sine- Gordon in ip2 which describes the commensurate-incommensurate 
transition. When 8f is small, the phase difference ip2 is zero an( i kinks are 
not favorable. This implies that there is no vortices (q = 1) or the vortex 
lattice is commensurate with the ladder (q > 2). Above a critical value, 
fife ~ \A^ 5 kinks appear separating commensurate domains and the vor- 
tex lattice is incommensurate (q > 2). For q = 1 this corresponds to flux 
penetration in the ladder [16]. The commensurate region decreases with 
increasing charging energy E c and vanishes at a critical value. Including 
vortices in 4>\ , fa will result in additional disordered phases with exponen- 
tially decaying correlation functions where the ladder has an insulating 
behavior. For different values of q the topology of the phase diagram is 
indicated in Fig. 2. 

The topology of the phase diagrams in Fig. 2 is strongly dependent 
on q for f = p/q, displaying direct vortex commensurate-incommensurate 
transitions in addition to superconductor-insulator transitions. Allowing for 
finite fugacities of the vortices may shrink to zero the intermediate phase at 
8f = for q = 3, as renormalization-group analysis suggests [11]. However, 
this should persist for q > 4. 

The intermediate incommensurate vortex phase at 8f = for q > 4 
should have interesting transport properties which should be in principle 
experimentally measurable. In fact, the vortex lattice is effectively depinned 




Figure 2. Qualitative zero-temperature phase diagram as a function of charging energy 
E c and small deviations 8f = f — f from a commensurate frustration f = p/q. C, IC 
and I denote commensurate-vortex, incommensurate-vortex, and insulating phases. For 
q = 1 vortices are absent in the C phase. 



by quantum fluctuations in this phase. If a driving current is applied trans- 
verse to the ladder, we expect that it will lead to vortex motion along the 
ladder and give rise to a finite resistance even at zero temperature. 

3. Critical behavior at / = 1/2 

For increasing charging energy the one- dimensional ladder undergoes a su- 
perconductor to insulator transition at zero temperature. As discussed in 
Sec. 2, the details of this transition will depend both on the value of the 
frustration parameter / and the deviation from commensurability 8f. As 
in the case of two-dimensional superconducting films [13] and Josephson- 
junction arrays[4, 2, 3], the universality class of this transition is a problem 
of great interest. 

For a chain of Josephson junctions the superconductor-insulator tran- 
sition has been identified with that of the KT transition[22]. Also, for the 
case of two coupled chains forming a ladder, in the absence of magnetic 
field, this critical behavior remains unchanged. This can been seem from 
Eq. (3) by taking q = 1 and 8f = 0, in which case the resulting Hamilto- 
nian has a KT transition when vortices are included [11]. However, in the 
presence of the magnetic field the situation is more complicated. 

For / = 1/2 and 8f = 0, the effective Hamiltonian of Eq. (3) is expected 
to have a critical behavior in the universality class of the two-dimensional 
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classical XY-Ising model [15, 18, 19] defined by the classical Hamiltonian 
(3H = - + a r a r t cos(0 r - 9 r ,) + Ca r a r ,} (4) 

rr' 

The existence of both XY and Ising-like excitations in this case result from 
the frustration and are associated with the continuous U(l) symmetry of 
the phases of the superconducting order parameter and the plaquette chi- 
ralities which measures the direction of circulating currents in the ladder. 
This model has a phase diagram consisting of three branches in the ferro- 
magnetic region. One of them corresponds to single transitions with simul- 
taneous loss of XY and Ising order. Further away from the branch point, 
this line of single transitions becomes first order. The other two lines cor- 
respond to separate KT and Ising transitions. The ladder at / = 1/2 is 
represented by a particular path trough this phase diagram which will dep- 
dend on the ratio E x /E y , between the interchain and intrachain Josephson 
couplings. As a consequence of this identification, the universality class of 
the superconductor-insulator transition should depend on the ratio between 
the coupling constants. 

The critical behavior of the ladder at / = 1/2 has been studied re- 
cently using numerical techniques[15, 19]. The ladder at this value of / 
can be regarded as a one- dimensional frustrated quantum XY (ID FQXY) 
model by choosing a particular gauge, A = xBy, for the vector potential 
in Eq. (1). In this gauge, the ladder of Fig. 1 is equivalent to an XY model 
with interchain coupling E x and intrachain couplings E y for one chain and 
— E y for the other. The critical behavior was studied using a path-integral 
Monte Carlo transfer matrix and finite-size scaling. In this formulation the 
ID FQXY model is mapped into a 2D classical statistical mechanics prob- 
lem where the ground state energy of the quantum model of finite length 
L corresponds to the reduced free energy per unit length of the classical 
model defined on an infinite strip along the imaginary time direction. The 
parameter a = ^/E y /E c plays the role of an inverse temperature in the 
2D classical model. The scaling behavior of the energy gap for kink exci- 
tations (chiral domain walls) of the ID FQXY model corresponds to the 
interface free energy of an infinite strip in the classical model. For large a 
(small charging energy E c ), there is a gap for creation of kinks in the anti- 
ferromagnetic pattern of chiralities Xp an d the ground state has long-range 
chiral order. At some critical value of a, chiral order is destroyed by quan- 
tum fluctuations with an energy gap vanishing as \a — a c \ v . At the critical 
point, the correlation function decay as power law with an exponent r]. The 
free energy per unit length /(a) on the infinite strip can be obtained from 
the largest eigenvalue A of the transfer matrix between different time slices 
as / = — In A . To obtain A , the Monte Carlo transfer-matrix method [20] 
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was used. Similar approach has been used to study the critical behavior 
of the 2D frustrated classical XY model [12]. The results for the critical 
coupling a c and critical exponents v and 7/ for two different values of the 
ratio E x /E y are indicated in Table 1. 

TABLE 1. Critical cou- 

plings (a c = (Ey/Ec) 1 / 2 ) and exponents 
(y, rj), obtained from finite-size scaling 
analysis of interfacial free energies for 
two values of the ratio between inter- 
chain and intrachain couplings (Ex/Ey) 



E X j Ey 




V 


n 


1 


f. 04(f) 


0.81(4) 


0.47(4) 


3 


1.16(2) 


1.05(6) 


0.27(3) 



For equal couplings E x = E y , the results for the critical couplings r/ 
and v differ significantly from the pure 2D Ising-model exponents (y = 
1,7/ = 0.25). In the XY-Ising model this correspond to a path through 
the phase diagram in the single transition region [18]. It is interesting to 
note that this result is also consistent with similar calculations for the 
2D frustrated classical XY model [12]. In general, the critical behavior of 
a d dimensional quantum model is in the same universality class of the 
d+1 dimensional classical version. This would suggest that the ID FQXY 
model is the Hamiltonian limit of the 2D classical model. Apparently, this 
is not the case although their critical behavior appears to be in the same 
universality class. 

For E x /E y = 3, the results indicated in Table 1 are in good agreement 
with pure 2D Ising values. In the phase diagram of the XY-Ising model, 
this correspond to a path in the decoupled region where KT and Ising 
transitions can take place at different points. In Fig. 3, the results for the 
helicity modulus are shown [19]. 

This quantity measures the response of the system to an imposed twist 
of 7r along the strip and is given by 7 = 2AF/ir 2 for large system sizes, 
where Af is the free-energy difference between strips with and without 
an additional phase mismatch. If the model were decoupled, the transition 
should be in the KT universality class, where one knows that the transition 
is associated with a universal jump of 2/ it in the helicity modulus [21]. The 
critical coupling can be estimated as the value of a at which AF = it which 
gives a c = 1.29. When compared with the critical coupling for destruction 
of chiral order in Table 1, a c = 1.16, it clearly indicates that there are two 
decoupled transitions. Somewhere in between a = 1 and 3 there should 




Figure 3. Interfacial free energy AF for a system of size L = 12 resulting from an 
imposed phase twist of ir. Vertical arrows indicate the location of the Ising and KT 
transitions and the horizontal arrow the value AF = ir from where the KT transition is 
located. 



be a bifurcation point where a single transition with simultaneous loss of 
phase coherence and chiral order decouples into separate transitions. Since, 
the superconductor to insulator transition is to be identified with the loss 
of phase coherence, in the decoupled region this transition is in the KT 
universality class but in the single region it is in a new universality class. It 
is interesting to note that, unlike the zero field case, the universality class 
of the superconductor-insulator transition at / = 1/2 depends on the ratio 
between the inter- and intra-chain couplings. 

4. Conclusions 

The zero-temperature phase diagram of a periodic Josephson-junction lad- 
der was considered. As a function of the magnetic field, the topology 
of the phase diagram is strongly dependent on q for f = p/q, display- 
ing direct commensurate-incommensurate vortex transitions in addition to 
superconductor-insulator transitions. For q > 4, an intermediate incommen- 
surate phase is also possible where the vortex lattice is effectively depinned 
by quantum fluctuations. Although these are zero-temperature properties, 
they can also be observed at finite but low temperatures when the thermal 
correlation length exceeds the sample size. 

The superconductor-insulator transition at / = 1/2 is found to depend 
on the ratio between interchain and intrachain coupling constants. The 
result is consistent with the predictions based on the phase diagram of the 
classical XY-Ising model. For increasing interchain couplings the XY and 



9 



Ising-like excitations decouple, giving rise to pure Ising critical behavior for 
the chirality order parameter and a superconductor-insulator transition in 
the universality class of the KT transition. For equal couplings, there is a 
single transition with simultaneous loss of phase coherence and chiral order 
with critical exponents in agreement with the finite-temperature transition 
of the frustrated XY model. 
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